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A method for realizing asymmetric (one-way) transmission of discretized light in modulated, linear and purely 
passive optical lattices is suggested, which exploits the idea of unidirectional coherent perfect absorption. The 
system consists of a linear photonic lattice of coupled resonators or waveguides, side coupled to a chain of lossy 
elements, in which light can avoid the occupation of the dissipative sites when propagating in one way, but not 
in the opposite one. Non-reciprocity requires modulation of the resonator/waveguide parameters, realizing a 
dissipative optical Aharonov-Bohm diode with non-reciprocal behavior. © 2015 Optical Society of America 
OCIS codes: 230.3240, 230.3120, 230.7370, 000.1600 


Unidirectional transport is at the heart of many funda¬ 
mental problems and applications in science and technol¬ 
ogy. In optics, unidirectional transport enables to realize 
an optical diode, which is generally obtained by time- 
reversal symmetry breaking using magnetic fields in lin¬ 
ear media or by exploiting nonlinear interactions. In the 
past few years, several efforts have been devoted toward 
the realization of miniaturized on-chip optical isolators 
using novel concepts and design methods. These include 
the use of integrated magneto-optics materials [ 1 , 2 ], 
time-varying optical structures [3-7], and 'PT-symmetric 
nonlinear optical media [8-11]. Asymmetric propagation 
has been also shown to occur for discretized light in pho¬ 
tonic lattices, i.e. arrays of coupled optical resonators or 
waveguides, where transport is described by an effective 
discrete Schrodinger equation [11-18]. In such studies, 
one-way transmission was realized by the introduction 
of some nonlinear and asymmetric defects in the lat¬ 
tice [11-13,16-18] or using magnetic media [14,15]. 



Fig. 1. (Color online) (a) Schematic of a linear array of 
coupled resonators or waveguides (C) side-coupled to a chain 
of (5 + 1) dissipative sites (A). The resonator frequencies are 
modulated as discussed in the text, (b) Effective static lattice 
model obtained in the rotating-wave approximation, with the 
appearance of a Peierls phase 0, corresponding to an effective 
magnetic flux <f> threading each triangular plaquette (see the 
inset), (c) Realization of the Peierls phase using an array of 
helical waveguides. 


In this Letter a different route toward the realization 
of asymmetric transport in discrete optical systems is 
proposed, which is based on the idea of unidirectional 
coherent perfect absorption (CPA). Unidirectional trans¬ 
parency is obtained by realizing a dissipative Aharonov- 
Bohm diode, in which an effective magnetic field is in¬ 
troduced by periodic modulation of the optical parame¬ 
ters [6,19-21]. The method works for linear and passive 
lattices and does not require nonlinear, active nor mag¬ 
netic media. 

We consider light propagation in a one-dimensional 
lattice C, which is side coupled to a number (5 + 1) 
of lossy (dissipative) sites A as schematically shown 
in Fig. 1 (a). Such a structure can be realized by cou¬ 
pled resonator optical waveguides (CROW) or by arrays 
of evanescently-coupled optical waveguides. Light trans¬ 
port is described by the coupled-mode equations 

dc S 

^~dt ~ UJn M Cn A ^( c n+1 + c n- 1) + & (1) 

1=0 

i~dt ~ — ^7 ) a l + o(ci + c/_|_i) + 0(ai+i + a/_i) (2) 

(1 = 0 , 1 ,..., 5), with ai = 0 for l < 0 and l > S. 

In Eqs.(l,2), uj n ,c and care the resonator frequen¬ 
cies (propagation constants) of the resonators (waveg¬ 
uides) in the main lattice C and in the side-coupled chain 
A, respectively, c n (t) and a n (t) are the respective field 
amplitudes, k, 9 and a are the coupling constants [as 
shown in Fig. 1(a)], and 7 is the loss rate in the lossy 
resonators (waveguides) A. Note that in CROW struc¬ 
tures t is the time variable, whereas in waveguide arrays 
it represents the longitudinal spatial propagation coor¬ 
dinate. For 5 = 0 and for uniform resonator frequencies 
^ n ,A = ^i,b = the structure of Fig. 1(a) is atypi¬ 

cal configuration that displays a Fano resonance in trans¬ 
mittance for uja 7 ^ However, for static resonators, i.e. 
ujc,a independent of time, the transmittance (as well as 
reflectance) is independent of the incidence side. Such 
a symmetric (reciprocal) behavior holds regardless the 
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number (S' + 1) of resonators in the chain A nor the 
value of the loss rate 7. To obtain a non-reciprocal be¬ 
havior, we introduce an effective magnetic field for pho¬ 
tons by proper modulation of the resonator frequencies 
uo n ,c and oj^a so as to introduce Peierls phases in the 
hopping rates [6,19-21]. Previous works have suggested 
and discussed in details different methods to realize an 
effective gauge field for photons in CROW and waveg¬ 
uide lattices [19-24]. Here we consider the method of 
a modulated gradient index ramp, which is inspired by 
the method of photon-assisted tunneling used to intro¬ 
duce synthetic gauge fields for matter wave systems [25]. 
Specifically, let us assume 


Transmittance Reflectance 



u n ,c(t) = u c +2nF x (t) , ui iA (t) = u A +(2l+l)F x (t)+F y (t) 

(3) 

where ujc , are the static resonator frequencies, 

F x (t) = —T x uj sm(u)t + cp) , F y (t) = TyUj cos(u:t + ip) 

( 4 ) 

are the sinusoidally-modulated gradient index in the hor¬ 
izontal (x) and vertical ( y ) directions, uj is the modula¬ 
tion frequency, and 

= r sm(<p) , Ty = r cos(<p) (5) 

are the modulation amplitudes. In the high-frequency 
modulation limit uo ft, cr, 0 and assuming the reso¬ 
nance condition Muj = ujc — u A for some integer M, 
an effective stationary lattice with a fictitious mag- 



wave number q wave number q 


Fig. 2. (Color online) Dispersion curves E±(q) (real and 
imaginary parts) of the two minibands for the binary lattice 
of Fig. 1(b) in the S oo limit. Parameter values are: fti = 
ft 3 = 1, ft 2 = 0.3, (j) = 7t/2, and 7 — 0.2. 


Fig. 3. (Color online) Spectral transmittance Tl,r(c[) = 
\tL,R(q )\ 2 and reflectance R(q) = |rz,(g)| 2 = |rft(g )| 2 for left 
(L) and right (R) incidence sides in the lattice of Fig. 1(b) 
with S = 0 and for a few values of ft 2 /fti: (a) ft 2 /fti = 1, (b) 
ft 2 /fti = 0.5 and (c) ft 2 /fti = 0.2. In the plots the loss rate 7 
is chosen according to the CPA condition (14). 

and where J n is the Bessel function of the first kind and 
of order n. Equations (6) and (7) describe the effective 
stationary lattice shown in Fig. 1(b), in which the Peierls 
phase mimics an effective magnetic flux threading each 
closed triangular loop. Such an effective magnetic flux, 
together with the dissipation in the lattice sites A, can 
provide asymmetric transmission along the main lattice 
C. To this aim, let us first consider the limiting case 
S —> 00, i.e. of two infinitely-extended side-coupled ar¬ 
rays C and A forming a binary lattice. The band struc¬ 
ture E = E(q) of the binary array can be readily ob¬ 
tained from Eqs.(6) and (7) with the Ansatz C n (t ) = 
Cexp[— iE(q)t — iqn \, A n (t) = Aexp[—iE(q)t — iqn \, 
where —7r < q < it is the Bloch wave number. The dis¬ 
persion relations for the two lattice minibands read 

E±(q) = (fti + ft 3 ) cosq - i- ( 10 ) 

± ]J [(fti - ft 3 ) COS q + + 2ft|[l + cos (q + <j>)\. 


netic flux is obtained (see e.g. [19, 20]). In fact, af¬ 
ter setting c n (t) = C n (t) exp [—2inQ x — ico^t], ai(t) = 
Ai(t) exp[—(2/ + 1 )iQ x — iO y — iw A t \, where & x (t) = 
Ta, cos (uit + cp) and Q y (t) = T y sm(u)t + cp), under the 
rotating-wave approximation Eqs.(l) and (2) reduce to 
the following coupled-mode equations for the amplitudes 
C n and Ai 

dC n ^ ^ 

i ~j7 L = K i( c n +1 + C n - 1 ) + ft 2 ^2(S n ,i exp+ 5 nA +i)Ai (6) 

1=0 

dAi 

i—— = —i'yAi + ft 3 (A + i + Ai _ x ) + ft 2 exp+ ft 2 ^+i (7) 
at 

(,l = 0,1, 2,.., 5), where we have set 

k 1 = ftJo(2rsin(^) , K 2 = crJ M <X) (8) 

ft3 = 6Jo(2Tsm(p) , cf) = 2M(p (9) 


A typical behavior of the real and imaginary parts of 
E±(q) are shown in Fig.2. Note that, for a small coupling 
ft2 between the two sublattices C and A, the miniband 
E+(q) corresponds to light localized mainly in the lat¬ 
tice C, whereas the miniband E-(q) corresponds to light 
being mostly localized in the (dissipative) lattice A. For 
4> 7^ 0 the symmetry of the dispersion relations under the 
change q -A — q is broken, which means that Bloch waves 
propagating in the forward (0 < q < tt) and backward 
(—7r < q < 0) directions in the lattice behave differ¬ 
ently. Interestingly, at q = qo = tt — (f) > 0 the forward 
propagating wave for the upper miniband (E + ) is not 
attenuated, i.e. lm(E+(qo)) = 0, whereas the backward 
propagating wave at q = —qo it is, i.e. unidirectional 
absorption is realized. Such a result indicates that the 
system behaves non-reciprocally for forward and back¬ 
ward wave propagation, and can be thus exploited to 
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Fig. 4. (Color online) (a-c) Behavior of the transmittance 
Tl,r (for left and right incidence sides) and reflectance R 
versus the loss rate 7 for increasing values of the number 
(5 + 1) of lossy sites in the chain C: (a) S = 2, (b) S = 4, 
and (c) S = 9. Parameter values are: ki = K 3 = 1, ^2 = 0 . 2 , 
and 0 = 7 r/ 2 . In (d) the behavior of the transmittance Tr(^) 
versus 7 is shown at q — 7t/2. 


realize an optical diode. 

Let us now consider the case of a finite number 
(S + 1) of dissipative resonators in the chain A. With¬ 
out loss of generality, let us assume > 0. The 
transmission (t) and reflection (r) coefficients for left 
(L) and right (R) incidence sides in the array C can 
be calculated by standard methods. For left-side inci¬ 
dence, the scattering solution is of the form C n (t ) = 
exp (—iqn — iEt) + 7 +(g, (j>) exp(iqn — iEt) for n < 0 
and C n (t) = tL(q, (f>) exp [—iq(n — S — 1 ) — iEt] exp(ign) 
for n > S + 1 , whereas for right-side incidence one has 
C n (t) = (f>) exp(ign — iEt ) for n < 0 and C n (t ) = 

exp[iq(n—S— 1)— iEt]-\-rR(q, (j>) exp[— iq(n — S— 1)— iEt] 
for n > 5 + 1. In such previous relations, we have set 
E = 2fti cos q and assumed 0 < q < n. One can show 
that the following general relations hold 

TR{qA) = r L {q,4>) , t R (q,<j>) =t L (q,-<l>). ( 11 ) 

This means that, while the system behaves sym¬ 
metrically in reflection, it behaves asymmetrically in 
transmission. A simple analytical form of the reflec¬ 
tion/transmission coefficients can be given in case of a 
single dissipative site. For 5 = 0 one has: 




rL{qA) = 


2 i sin q [/ci (E + ij) + /c| exp(i^)] 

/ci (E + Z7)[exp(2zg) — 1] — 2/c| cos 0 — 2/c| exp (iq) 
2 /c| (cos q + cos <j>) 

ki (E + Z7)[exp(2zg) — 1] — 2/c| cos</> — 2/c| exp(zg) 



0 

lattice site n 


lattice site n 


Fig. 5. (Color online) Propagation of a Gaussian wave packet 
of carrier wave number q — 7r/2 (snapshot of |c n (t )| 2 in a 
pseudocolor map), for left and right incidence sides, in the 
lattice C of Fig. 1 (a) with modulated parameters and with 
5 = 0, for (a) lj/k = 8 , and (b) uj/k — 20. The other pa¬ 
rameter values are given in the text. In (a) the estimated 
transmittances are Tl — 89.4%, Tr ~ 16%, whereas in (b) 
one has Tl — 95%, Tr ~ 4.8%. 


where E = cos q. In the following we will consider 
the case 0 = 7t/2, which is the most interesting one, and 
will omit the dependence of the reflection/transmission 
coefficients on <fi. From Eqs.(12) and (13) one has |£l (<7 = 
7r/2)| = 1, | t R {q = 7r/2)I = |7«i - K 2l/(7«i + K l) and 
\^l,r{q = 7 r/ 2 ) | = 0. In particular, provided that the loss 
rate 7 satisfies the condition 


7 = «i/ K i > 


(14) 


( 12 ) 


(13) 


unidirectional CPA is obtained: namely, a light wave 
propagating in the forward direction is fully transmitted, 
without being reflected nor absorbed, whereas a light 
wave traveling in the backward direction is neither re¬ 
flected nor transmitted, rather it is fully absorbed in the 
site A. Hence, when the condition (14) is met, the optical 
system realizes perfect optical isolation at q = 7r/2 be¬ 
cause of unidirectional CPA. We not that, while CPA in 
linear dissipative optical structures generally shows a re¬ 
ciprocal (symmetric) behavior [26-28], i.e. it is indepen¬ 
dent of the incidence side, in our lattice system CPA is 
a unidirectional phenomenon because of the presence of 
the Peierls phase 0. The structure of Fig. 1(b) with 5 = 0 
basically realizes a discrete and dissipative Aharonov- 
Bohm interferometer for light waves, allowing for the 
nonreciprocal transmission of the structure. A typical 
behavior of the transmittance Tji^(q) = \tL,R.(q )\ 2 and 
reflectance R(q) = \TR^{q )\ 2 for 5 = 0 and for a few 
values of ^ 2 /^ 1 , corresponding to the CPA condition 
(14), is shown in Fig.3. Note that, as ^ 2/^1 decreases, 
the bandwidth around q = 7 r /2 of unidirectional CPA 
shrinks. Hence, to realize a relatively broadband optical 
isolation, a strong coupling of defect A with the main 
lattice C is required. 

For an arbitrary number (5 + 1) of dissipative sites in 
the chain A, the condition for unidirectional CPA devi¬ 
ates from Eq.(14), and one has to resort to a numerical 
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analysis to determine the optimum value of 7 . However, 
as S increases, CPA is observed for a wide range of 7 
values. As an example, Fig.4 shows the behavior of trans¬ 
mittance and reflectance versus the loss rate 7 for a few 
values of S. Interestingly, even for a small value of ^ 2 /^ 1 , 
broadband unidirectional CPA can be realized increasing 
the value of S [see Fig.4(c)]. 

Non-reciprocal transmission predicted by the effec¬ 
tive lattice model of Fig. 1(b), i.e. in the rotating-wave 
approximation, has been confirmed by beam propaga¬ 
tion simulations within the original modulated lattice 
model [Eqs.(l) and (2)]. As an example, Fig.5 shows the 
numerically-computed propagation of a Gaussian wave 
packet, for either left or right incidence sides, with carrier 
wave number q = it/2 and for parameter values 5 = 0 , 
a /k = 1.6, M = 1, 9 ? = 7r/4 (corresponding to the mag¬ 
netic flux <j) = 7r/2), r = 1, w/k = 8 in Fig.5(a), and 
uo/k = 20 in Fig.5(b). The loss rate y/zs = 0.88 has 
been tuned to satisfy the CPA condition (14), where the 
values of k\/k ~ 0.56 and — 0.7 are estimated 

using Eq.(8). The simulations clearly show strong asym¬ 
metric transmission, even for a relatively low modulation 
frequency [Fig.5(a)]. 

Finally, let us briefly discuss a possible experimental 
realization of asymmetric transmission exploiting unidi¬ 
rectional CPA and based on spatial light transport in 
lattices of dielectric optical waveguides [29]. In such an 
optical system, the Peierls phase can be implemented 
by periodically rotating the optical axis of waveguides 
along the propagation direction [22,24,29]; see Fig.1(c). 
Taking an helical path of the waveguide axis with spatial 
frequency uj = 27r/A described by the parametric equa¬ 
tions xo(t) = X 0 sin (uit + ip) and yo(t) = —To cos( cj£+ (/?), 
the effective index gradients given in Eqs.(3) and (4) are 
realized with T x = 27rn s aXoUJ /A and T y = 2irn s bYoU)/\, 
where A is the optical wavelength (in vacuum) of propa¬ 
gating waves, a and b are the horizontal and vertical dis¬ 
placements of waveguides in the lattice [see Fig.1(c)], and 
n s is the substrate (cladding) refractive index [29,30]. 
For example, let us consider propagation of visible light 
(A = 633 nm) in dielectric helical waveguides manufac¬ 
tured in fused silica by femtosecond laser writing [24,30]. 
Assuming n s ~ 1.46, 2 a = 16 /im, b = 4 /mi, k ~ 
1.5 cm -1 , the simulation of Fig.5 (a) corresponds to a he¬ 
lix with period A 0.52 cm and amplitudes Xo ~ 5 //m, 
Tq ~ 10 /mi, a loss rate 7 ~ 1.33 cm -1 (i.e. ~ 11.5 
dB/cm), and a full propagation length t ~ 10 cm. The 
detuning ujc — ^a = cj of the modal propagation con¬ 
stant of waveguide A from the main array C corresponds 
to a refractive index change 5n ~ A/A ~ 1.21 x 10 -4 of 
the core of waveguide A, which can be obtained by con¬ 
trol of the laser writing speed in the sample. Finally, the 
required loss rate in waveguide A can be realized, for 
example, by proper waveguide segmentation. 

In conclusion, a method for realizing asymmetric (one¬ 
way) transmission of discretized light in modulated op¬ 
tical lattices has been suggested, which is based on the 
idea of unidirectional coherent perfect absorption. The 


method works for linear and passive lattices and does 
not require nonlinear, active nor magnetic media. 
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